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q ■ Abstract 

f— v | We show, that higher analogs of the Willmore functional, defined 

C"-- ' on the space of immersions M 2 — > M 3 , where M 2 is a two-dimensional 

torus, M 3 is the 3-dimensional Euclidean space are invariant under 
conformal transformations of K 3 . This hypothesis was formulated re- 
cently by I. A. Taimanov. 
Pl ^ 1 Higher analogs of the Willmore functional are defined in terms of 

the Modified Novikov-Veselov hierarchy. This soliton hierarchy is asso- 
ciated with the zero-energy scattering problem for the two-dimensional 
/\i ■ Dirac operator. 

1 Introduction 

To start with, we would like to recall the following interesting fact from 
the theory of 2-dimensional surfaces in R 3 (see ]2(|, p. 110 and references 



&jd: 



therein). Let X : M 2 — > M 3 be a smooth immersion of a compact orientable 
surface M 2 into the Euclidean space IR 3 (i.e. a smooth map from M 2 to R 3 
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such that its Jacobi matrix is non-degenerate everywhere on M 2 , but the 
image is allowed to self-intersect). Let T be the Willmore functional 

T= [ H 2 dS, (1) 

Jm 2 

where H denotes the mean curvature, dS is the volume element on M 2 gener- 
ated by the immersion. Then T is invariant under conformal transformations 
of R 3 . 

It is natural to pose the problem of constructing other conformal invariant 
functional of immersions. 

It is well-known, that many constructions from the soliton theory have 
natural analogs in geometry and vise versa. In particular, important informa- 
tion how to study immersions of 2-dimensional surfaces into R 3 using soliton 
methods can be found in the article |IJ by A. I. Bobenko. Any immersed 
surface possesses (at least locally) a conformal coordinate system (see for 
example M p. 110), i.e. a coordinate system such that ds 2 = f(z,z)dzdz. 
In conformal coordinates this immersion can be locally represented by the 
Generalized Weierstrass Formulas (see Section || below) and the potential 
U(z, z) is uniquely defined. In |T7| it was shown, that any analytic immer- 
sion of a compact orientable 2-dimensional manifold into R 3 can be globally 
represented by these formulas. 

The Generalized Weierstrass Formulas are based on the zero-energy eigen- 
functions of the two-dimensional Dirac operator with a real potential U(z, z). 
The zero-energy spectral problem for this operator arose in the soliton theory 
as an auxiliary linear problem for the hierarchy of Modified Novikov-Veselov 
Equations (MNV) (see M). These nonlinear integrable equations with 2 
spatial variables, introduced by L. V. Bogdanov ,have infinitely many con- 
servation laws. 

In B. G. Konopelchenko and I. A. Taimanov showed that the quadratic 
MNV conservation law 

#1 = 4/ / U 2 (z, z)dx Ady, z = x + iy (2) 



coincides with the Willmore functional T. In ]T8) I. A. Taimanov formulated 
a hypothesis that all higher MNV conservations laws also generate functionals 
on immersions of closed orientable 2-dimensional surfaces into R 3 , invariant 
under conformal transformations of R 3 . He also did some numerical experi- 
ments with surfaces of revolution, confirming this assumption. 



During his visit to the Freie Universitat, Berlin in September-October 
1996, I. A. Taimanov attracted the authors attention to this problem. In the 
present text we prove this hypothesis for immersions of tori into M. 3 . 

Of course, it would be natural to extend this result to immersions of 
higher genus surfaces into R 3 . But here we meet the following problem. In 
contrast with Hi the higher conservation laws are non-local in terms of the 
potential. For double-periodic potentials they are defined in terms of the 
zero-energy dispersion curve (the Riemann surface of the zero-energy Bloch 
function). If the genus is greater than 1, we have to study modular invariant 
Dirac operators in the Lobachevskian plane with a non-abelian group of 
translations. The corresponding Bloch theory has not been constructed until 
now, thus we could not define the higher conservation laws. It would be 
interesting to develop the corresponding Bloch theory, but this problem looks 
rather non-trivial. Thus we restrict ourself to the genus 1 case only, where 
all integrals of motion are well-defined. 

The idea of our proof is the following. We show, that infinitesimal confor- 
mal transformations of R 3 correspond to infinitesimal Darboux transforma- 
tions of the Dirac operator (infinitesimal dressings with degenerate kernels). 
It is convenient to express these deformations in terms of Cauchy-Baker- 
Akhiezer kernels, introduced by A. Yu. Orlov and one of the authors in || 
(see also the review ). From the explicit formulas for such deformations 
it follows, that the deformed zero-energy Bloch function is meromorphic on 
the same Riemann surface as the original one. Recalling, that this Riemann 
surface completely determines all conservation laws, we complete the proof. 

Remark 1 An alternative proof of the theorem, namely that conformal trans- 
formations do not change the zero-energy Bloch variety, was obtained by 
U. Pinkall (private communication). He calculated the action of finite con- 
formal transformations on the Bloch function using a quaternionic represen- 
tation of the Generalized Weierstrass formulas as suggested by G. Kamberov, 
F. Pedit and U. Pinkall in ||/. 

Remark 2 Finite Darboux transformations for the Dirac operator are dis- 
cussed in the book by V. B. Matveev and M. A. Salle |73|/ , Laplace transfor- 
mations for the Dirac operator are discussed in the paper by E. V. Ferapontov 

The authors are grateful to I. A. Taimanov for interesting discussions. 



2 Generalized Weierstrass construction 

Let L be the two-dimensional Dirac operator 

L- 



d z -U(z,z) 
U(z, z) d- z 



(3) 



with a real potential U(z,z). Let ty(z,z) be a zero-energy solution of the 
Dirac equation 



L^(z,z) = 0, V{z,z) 



*i(z,z) 
*2(M) 



(4) 



Then the Generalized Weierstrass Formulas (see [I7J and refs. therein) 
Xi(z, z) + iX 2 (z, z) = C 1 + tC 2 + if (mz', z')dz' - f 2 (z', z')dz') 



X x {z, z) - iX 2 (z, z) = C 1 -iC 2 + iJ {*\{z!, z')dz' - #?(*', z')dz') (5) 



Z|> 



X 3 {z, z) = C 3 -j (Mz', z')Vi{z', z')dz' + V 1 ( Z ', z')^^', z')^' 

defines a map of the plane K 2 to the Euclidean space R 3 . In ([5]) z is a fixed 
point in the z-plane and the integrals are taken over some path connecting 
the points z and z. From (|j) it follows, that the integrands in fl5|) are closed 
forms, thus the map does not depend on a specific choice of the path. Here 
Ci, C 2 , C-s are arbitrary real constants. 
The formulas (HD are equivalent to: 



d [a 2 X x + a x X 2 - a 3 X 3 ] 



y x (z,z) V 2 (z,z) 
-V 2 {z,z) Vi(z 7 z) 



dz 
dz 



#i(M) -^ 2 (z,z) 
V 2 {z,z) $i{z,z) 



(6) 



where a , er , a are the standard Dirac matrices 



(7 
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The Generalized Weierstrass Map is conformal, i.e. the metric ds 2 on R 2 
induced by this map is proportional to the standard one: ds 2 = g(z,z)dzdz. 

Assume that we have a map of a 2-dimensional torus into R 3 . Then the 
corresponding potentials U(z, z) is periodic 

U{z + f 1 ,z + T x ) = U{z + T 2 ,z + T 2 ) = U{z, z). (8) 

Also the eigenfunction y£?(z, z) is periodic or anti-periodic, i.e. 

^(z+T 1 ,z+f 1 ) = Wi${z,z), V(z+T 2 ,z+f 2 ) = m$(z,z), Wl = W 2 2 = 1. 

(9) 
The coordinate z is defined up to linear transformations z — > az + b, 
a, b G C, a ^ 0. Thus without loss of generality we may assume 

T x = 1, T 2 = r, Imr > 0. (10) 

Conditions (|8|) and (Q) are necessary, but, of course, not sufficient for 
periodicity of the Generalized Weierstrass map. Necessary and sufficient 
conditions for periodicity can be formulated in terms of the Bloch variety. 
They are obtained in a forthcoming paper by I. A. Taimanov and one of the 
authors (M.S.). We do not use these conditions in our text, thus we will not 
discuss them in further details 

3 Bloch function and Bloch variety. 

In this Section we assume that U(z, z) is real, smooth, and double-periodic 
(H). With any such potential we associate a one-dimensional subvariety Y in 
the two-dimensional complex space (C\0) 2 . 

The first object we need is the Bloch function. By definition, the Bloch 
functions if;(wi,W2, z, z) are quasiperiodic solutions of the Dirac equation ([|) 
with the following periodicity properties: 

l/j[w 1 ,W 2 , Z + 1,Z+1) = W 1 lp[w 1 ,W 2 , Z, z), , n > 

tp(Wi, W 2 ,Z + T,Z + f) = W 2 1p{w 1 ,W 2 , Z,z). 

The pairs of multipliers w±, w 2 possessing at least one non-zero Bloch 
solution form a complex one- dimensional subvariety T e (C\0) x (C\0) (see 



[|T^j). This variety is called the Bloch variety or the zero-energy disper- 
sion curve. For a generic potential U(z, z) the genus of T is infinite. 

The Bloch functions form a one-dimensional holomorphic bundle over T 
(it is shown below, that for generic A G T a Bloch solution is unique up to 
normalization). It is convenient to fix a section of this bundle ip(\,z,z), by 
assuming 

iP 1 {\,z,z)+i> 2 {\,z,z)\ z=Zl = l, (12) 

where z\ is an arbitrary fixed point. 

The logarithms of the multipliers Wi(A), u> 2 (A) 

p 1 (A) = -lnu; 1 (A), p 2 (A) = J-lnw 2 (A), (13) 

% i \t\ 

are called quasimomentum functions. Of course, they have non-trivial 
increments while going along cycles in T, and they are defined up to adding 
27mi, 2im 2 / \t\ respectively, where ri\ and n 2 are some integers. Thus the 
functions Impi(X), lmp 2 (\) are single- valued in T. The differentials of the 
quasimomentum functions 

dpi(X) = -^pi(\)d\, dp 2 (X) = ^P2(\)d\, (14) 

are single- valued and holomorphic on the finite part of T. 

In our text the Dirac operator @ is symmetric and the potential U(z, z) 
is real. Let us show, that the corresponding Bloch variety V possesses 7L 2 x 7L 2 
as a group of symmetries. An analogous statement for the fixed-energy Bloch 
variety corresponding to a two-dimensional self-adjoint Schrodinger operator 
was proved in [Tl| . The proof from [|TUJ may be applied to (|3[) after a minimal 



modification. 

The operator (||) with real potential U(z, z) has the following symmetry. 
If if;(wi,W2, z, z) is a Bloch solution of ([|), then the function 

V K ' ' ' J V -Mwi,w 2 ,z,z) ) y ! 

is also a Bloch solution of (||) with multipliers w\ and w 2 respectively. Thus 
the surface V possesses an antiholomorphic involution (we denote it by <j9 
for historical reasons) 

aO : T — ► T, o6 : (wi,w 2 ) — > (wi,w 2 ), (16) 



and 

ft(\,z,z)=n-\\)i(ad(\),z,z), (17) 

where n(\) is a scalar function, meromorphic in A and independent on z, z. 
It is less trivial to see that the surface T possesses a holomorphic involu- 
tion: 

a : T — > r, a : (wi,w 2 ) — > (lujf 1 ,^ 1 ). (18) 

To prove it, let us fix a generic point A G T. Denote by C WlW2 the Banach 
space of all locally square-integrable two-component complex-valued vector- 
functions on R 2 with the periodicity properties (|TTf) - The space £^-1^,-1 is 

naturally dual to C WlW2 . Namely, if f(z,z) G £^ ltU2 and g(z,z) G ^-i^-i, 
then we define a scalar product by 

1 1 

< /, g >=Sfdt 1 A dt 2 [/i(ti + rt 2 , ti + rt 2 )^i(ti + rt 2 , £1 + ft 2 )+ /-, ^ 

l iy J 

+/ 2 («i + rt 2 , t x + ft 2 )g 2 (tx + rt 2 , tx + rt 2 )] . 

Let f^ = ip(wx, w 2 , z, z), f^\ . . . , p n \ ... be the Jordan basis for the Dirac 
operator L in the space C WlW2 . Also let (f°\ if 1 ', . . . , <f n \ . . .be the dual 
basis in C w -i w -i. The functions <ft n > form a Jordan basis for the transposed 
operator L T . But L is symmetric with respect to this scalar product thus it 
has a zero eigenfunction in the space ^-i^-i. Hence if (wx,w 2 ) G T, then 
cr(wx,w 2 ) = (w^ 1 ,w 2 1 ) G T. 

One of the main properties of the Bloch variety T is the following: T may 
be treated as a complete set of integrals of motion for the Modified 
Novikov-Veselov hierarchy. 

Indeed, consider the space of all real- valued smooth double-periodic func- 
tions on C 1 = R 2 with a fixed pair of periods 1 and r. The Modified Novikov- 
Veselov hierarchy (MNV) (see Section [| below) defines an infinite collection 
of flows on this space 

dU(z,z,t 2n+1 ) 



dt 2n+ x 



K 2n+ x[U]+K 2n+1 [U], (20) 

^ ( y' W) = i (K 2n+l [U] - K 2n+l [U\) , (21) 

ot 2n+1 v ' 

where K 2n+ x[U] is some integro-differential operator in z, z. Here t 2n+ i, t 2n +i 
are parameters of these flows. 

7 



Statement 1 Let U(z,z,t2 n +i) be a solution of one of the MNV equations; 
let L(t2n+i) be the corresponding two-dimensional Dirac operator (j3[) depend- 
ing on an extra parameter t 2n +i ; let r(t 2 n+i) be the corresponding family of 
Bloch varieties. 

Then r(t 2n+1 ) = T does not depend on the MNV time £271+1 ■ 

Remark 3 Here and below we use the following notational convention. If 
we have a complete proof of a mathematical result we call it Theorem or 
Lemma. If we do not have a complete strict proof yet we use the word 
Statement. 

An analogous statement is well-known for soliton systems with one spa- 
tial variable. In Section [|we prove this fact at least for algebraic-geometrical 
potentials, corresponding to varieties T of finite genus. (Sometimes such po- 
tentials are called finite-gap potentials). It is rather clear, that our proof can 
be extended to all smooth potentials, but to do such extension strictly we 
need more detailed information about analytic behavior of the Bloch func- 
tions near infinity in the momentum space, than we have now. An appropri- 
ate analytic lemma for the one-dimensional Dirac operator, corresponding to 



the surfaces of revolution, was proved by one of the authors (M.S.) in [16 . 

There exists a different way (may be a more natural one) to get a strict 
proof of Statement |l|. It would be interesting to prove the following approx- 
imation property: 

Conjecture 1 Any smooth potential can be approximated by the algebraic- 
geometrical ones with the same periods. 

From such a result it would follow, that we can restrict ourself to the 
algebraic-geometrical potentials in our calculations. 

We have a map U(z,z) — > T[U] from the space of double-periodic real 
smooth functions to the space of complex subvarieties in (c\0) 2 , which is 
invariant under the whole MNV hierarchy. This map generates an infinite 
family of MNV "normal" conservation laws. Namely, let w\ be a generic 
point in C\0. Than we have an infinite collection of numbers w 2 [U] such 
that {w\i w 2 W]) £ r[[/]. From Statement p] it follows that these functionals 
w 2 [U] are conservation laws of the whole Modified Novikov-Veselov 
hierarchy. 



The functionals w 2 [U] are essentially nonlocal. In Section [| we show 
that under the same assumptions as in Statement [l] we can expand these 
functionals in some asymptotic series near infinity and the expansion coeffi- 
cients give us the standard "quasi-local" conservation laws. 

4 Conformal transformations of the Euclidean 
space R 3 , and MNV integrals of motion. 

In the previous Section we have associated with any double-periodic smooth 
real potential a Bloch variety T[U]. The map is constant on the trajectories 
of the MNV hierarchy. In this Section we associate with any immersion of a 
torus into R 3 a Bloch variety T and show, that F is invariant under conformal 
transformations of R 3 . 

Let M 2 be a torus with a fixed basis of cycles a, b. Let X : M 2 — > M 3 be 
a smooth immersion of M 2 into the Euclidean space. The standard metric 
on M 3 induces a conformal structure on M 2 . Let z be a conformal global 
coordinate on the universal covering space of M 2 . The coordinate z is defined 
uniquely up to affine transformations z — > cz + d. If we assume, that the 
shift of M 2 along the a-cycle corresponds to the shift z — > z + 1 then the 
coordinate z is defined uniquely up to shifts 

z^z + d. (22) 

The immersion X and coordinate z define a potential U(z, z) therefore 
also a Bloch variety T[U]. T[U] is invariant under the shifts ( p2]) thus it is 
completely determined by the immersion X and the cycles a, b, and we may 
write T[X, a, b]. 

Conformal transformations of the Euclidean space R 3 do not affect the 
conformal structure of M 2 , thus they leave the coordinate z invariant up to 
the shifts (0). Without loss of generality we shall assume that conformal 
transformations of R 3 do not change z. 

Now we are in position to formulate and prove our main result: 

Theorem 1 Let X : M 2 —> R 3 be an immersion of a torus with a fixed basis 
of cycles a, b into the Euclidean space; let T[X,a,b], be the corresponding 
Bloch variety. Then T[X,a,b] is invariant under conformal transformations 

ofm 3 . 

9 



Proof of the Theorem: 

Step 1: To start with let us recall the well-known facts about the group 
of conformal transformations of the standard Euclidean metric on R 3 (or on 
the sphere S 3 ) (see for example 0). This group is generated by the following 
transformations: 

1. Translations Xj — > Xi + (X )i. 

2. Rotations X -> AX, A e SO(3). 

3. Dilations X -> kX, k e K 

4. Inversions 

X % -> J^Z^— ■ (23) 

< -A — -<^0j -^ — ^-0 ^ 



5. Reflections 



X ^ X -2v<v,X >, <v,v>=l. (24) 



The connected component of the identity of this group is isomorphic to 
50(1,4) (see ||, page 143). The corresponding Lie algebra is generated by 
the following basis of infinitesimal transformations: 

1. Translations P a : 5Xi = 5{ a . 

2. Rotations fl a b, a < b: 5Xi = 5ibX a — S ia X b . 

3. Dilation D: SX t = X % 

4. Inversions K a 

3 

5X t = 2X { X a - S ia Y^ XjXj. (25) 

i=i 

Step 2: Let us calculate deformations of the potential U(z, z), the eigen- 
function *&(z, z) and the constants Cj in formulas ([5]) corresponding to all 
infinitesimal conformal transformations of R 3 . 

1. Translations simply shift the constants Cj and change neither U(z,z) 
nor ty(z,z). Thus they do not change T[X,a,b], and without loss of 
generality we can assume 

Cj = 0, j = 1,2,3. (26) 

10 



2. A simple direct calculation based on thr representation (^[) (we do not 
like to reproduce it here) shows that the rotations in R 3 correspond to 
the following transformations of the eigenfunction *£?(z, z) 



^ * 2 («, z) J \ V 2 {z, z) J^ H \-^ 1 (z,z) )' ] ' T IPI 

where a and /3 are some complex parameters. 
Both functions ^(z, z) and 



(27) 



* + <-M -!!(A)- (28) 

satisfy the Dirac equation (|4|) with the same potential U(z,z). Thus 
the rotations does not change the potential and T[X, a, b]. 

3. The dilation is generated by the scaling transform 

8$(z,z) = ±$(z,z) (29) 

and changes neither U(z, z) nor T[X, a, b]. 

4. The only nontrivial transformations of the Dirac operator correspond to 
the inversion generators. Up to conjugations by rotations all generators 
(|25| ) are equivalent. Thus it is sufficient to prove, that T[X, a, b] is 
invariant if we apply the generator —K 3 : 

5X 1 = -2X X X 3 , 5K 2 = -2X 2 X 3 , 6X 3 = -X 3 2 + X\ + X\. (30) 

Let us introduce the following notation 

W{z,z)=X 1 {z,z)-iX 2 {z,z). (31) 



The transformation ( p0[) corresponds to the following transformation 
of the function *£?(z, z) 

8V x (z, z) = -X 3 (z, z)#i(z, z) + iW(z, z)V 2 (z, z) . , 

6* 2 {z, z) = -X 3 (z, z)* 2 (z, z) - iW(z, z)V 1 (z, z) {6 > 

11 



Let us check it. 
SX 3 (z, z) = -j \(5^ 2 (z', z')^ 1 (z', z>) + V 2 (z', z')S^ 1 (z', z 1 )) dz' + 

Z(> 

+ (8V 1 (z', z')V 2 (z', z>) + ^ 1 (z', z')8$ 2 (z', z')) dz' = 

= - J U-2X 3 (z>, z')Mz', z')*i(z',z')- 

-iW(z', z')^l{z', z!) - iW(z', z')^l(z', z 1 )) dz'+ 

+ (-2X 3 (z>, z'^^z', z')V 2 (z', z')+ 
+iW(z', z')mz', z!) + iW(z', z')^j(z', z')) dz' = 

= -f{(2X 3 (z',z')(d z ,X 3 (z',z'))- 

-W{z', z')(d z ,W(z', z')) - W(z', z'){d z ,W{z', z')j) dz'+ 

+ (2X 3 (z',z')(d- z ,X 3 (z',z'))- 
-W(z', z')(d z ,W(z', z')) - W(z', z')(d z ,W(z', z'))) dz' 

= -J \(d z ,Xl(z',z') - d z ,(W(z',z')W(z',z'))) dz' + 
+ (d- zl Xl(z', z') - d z ,(W(z', z')W(z', z'))) dz' = 
= W(z,z)W(z,z)-X 2 3 (z,z) 
SW(z, z)=iJ2 \8V 2 (z', z')^ 2 (z', z')dz' - S^^z', z')®i(z', z')dz' 



J 2 ( - X 3 (z', z')W 2 2 (z', z') + W(z', *')*i(*'> z')^ 2 (z', z')\ dz' + 

20 L ^ ' 

+ (X 3 (z', z')W 2 (z', z') + W(z', z')$ 2 {z', z')^ 1 (z', z')) dz' = 
2J\(- X 3 (z', z')(d z ,W(z', z')) - W(z', z')(d z ,X 3 (z' : z'))) dz' + 

20 LV ' 

+ ( - X 3 (z',z')(d- z ,W(z',z')) - W(z',z')(d z/ X 3 (z',z'))) dz' = 

z 

(d zl (X 3 (z',z')W(z',z'))) dz' + ((d- z/ (X 3 (z',z')W(z',z'))) dz' 



20 



= -2W(z,z)X 3 (z,z) 
The corresponding transformation of the potential U(z, z) reads as 

SU(z, z) = V^z, z)®i(z, z) - V 2 (z, z)V 2 (z, z). (33) 



12 



Step 3: Let us calculate the the deformation of the Bloch functions 
corresponding to (|33|). 

Let ip(X, z, z) be the Bloch function of L. For any A such that at least 
one of the functions Imp x (\), Imp,, (A) is not equal to zero ("non-physical" 
A) define the following pair of functions 

fix(A, z,z) = J i/> 2 {\, z', z')$i(z', z')dz' + V>i(A, z', z')$ 2 (z', z')dz' 

f (34) 

fi 2 (A, z,z) = J V 2 (A, z', z')^ 2 (z', z')dz' - V> X (A, z', z')^x{z\ z')dz> '. 



where the integrals are taken along an arbitrary path in the z-plane, connect- 
ing the points z and oo such that the integrand decays exponentially along 
this path. The integrands in ( [51D are closed 1-forms thus the integrals do 
not depend on a concrete choice of the path. Using the same arguments as 
in Section |T2| below we may easily prove that f2i(A, z, z) and fi 2 (A, z, z) are 
meromorphic in A outside infinity. 

The function *&(z, z) is double-periodic or anti-periodic in z (see (^)), thus 
the functions Oi(A, z, z), S1 2 (A, z, z) have the following periodicity properties: 

Q k {X,z + l,z + l) = W 1 w 1 (X)Q k {X,z,z) 

Q k {X,z + T,z + f) = W 2 w 2 (X)Q k (X,z,z) /c " i ' / - V°) 

(see (|) for the definition of Wi, W 2 .) 

Lemma 1 The variation of the Bloch function ip(X,z,z) corresponding to 
flffl) reads as 



8ipi(X, z, z) = Oi(A, z, z)i£?i(z, z) - fi 2 (A, z, z)^ 2 (z, z) + a(A)V>i(A, z, z) 
5i/j 2 (X, z, z) = Oi(A, z, z)^ 2 {z, z) + S7 2 (A, z, z)\j>i(,2, z) + a(A)-^ 2 (A, z, z) 

(36) 
where a(A) zs some meromorphic function, fixed by the normalization condi- 
tion: 

5M\z,z) + 5M^^z)\ z=Zi =0, (37) 

Proof of Lemma [l]. To start with, let us recall a simple fact from Bloch 
theory (see for example [|Kj). 

Generically, if we calculate variations of the Bloch function, we deform 
T, and we can not assume both 5wi(X) = and Sw 2 (X) = simultaneously. 

13 



To compare functions on different subvarieties in C 2 we have to fix some 
connection. The simplest way to do this is to assume 6wi(X) = 0. 

Then the variation of the Bloch function can be found as the unique 
solution of the linearized Dirac equation 

8L$(\, z, z) + L6ip{\, z,z) = (38) 

satisfying fl3T|) such that 

6iJ(X,t 1 + Tt 2 ,h + ft 2 ) = 0([1+ \t 2 \] e ^i^+^^A . (39) 



A simple direct calculation shows that flgq ) solves flog). From fl35|) and 
@ it follows, that 

Sipjwx, w 2 ,z + l,z + l) = wiSt/Awi, w 2 , z, z), ,^ . 

8ip(wi, W 2 , Z + T, z + f) = w 2 Sip(wi, w 2 , z,z), 

thus variations of the type fl36|) satisfy (|39|) . This completes the proof. 

The function 5-0(A, z, z) defined by (|3~6l ) has the same periodicity prop- 
erties as the original Bloch function ip(X,z,z) (see formulas (|TTD and (^0|) 
respectively). Thus our special variations satisfy Swi(X) = and 5w 2 (A) = 
simultaneously. 

Step 4: From ( |36l) it follows, that if we apply infinitesimal conformal 
transformation — K 3 to our immersion, the Bloch functions of the deformed 
Dirac operator are meromorphic on the same variety T[X, a, b] as the original 
Bloch functions and have the same multipliers W\, w 2 . Thus our deformation 
does not change T[X, a,b}. This completes the proof of Theorem [l|. 

Example 1 Surfaces of revolution. Let 7 be a closed non-selfintersecting 
curve in the half-plane X 2 = 0, Xi > in K 3 . Rotating 7 about the axes 
Xi = X 2 = we get a surface of revolution. It is always a torus with a fixed 
pair of periods. 

Such surfaces are essentially simpler from the soliton point of view. Po- 
tentials U(z, z) corresponding to such surfaces depend only on one real vari- 
able x = Rez. Instead of the fixed energy spectral transform for the 2- 
dimensional Dirac operator we have the spectral transform for the 2x2 first- 
order matrix differential operator in one variable. Periodic direct spectral 
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transform for such operators (for both finite-gap and infinite-gap potentials) 
was developed by one of the authors (M.S.) in fiTb}] . 

The MNV equations for surfaces of revolution are reduced to the well- 
studied Modified Korteweg-de Vries equations (MKdV). In contrast with MNV 
all higher MKdV conservation laws are local in terms of the potential. Equa- 
tion t ffll ) in this situation was integrated by V. K. Melnikov in £Q{] in the 
class of potentials sufficiently fast decaying at infinity. Our theorem for the 
surfaces of revolution does not follow formally from / |7^|/ because the periodic 
MKdV theory and the decaying at infinity one use different technical tools. 
Nevertheless it is possible to simplify essentially our proof in this specific 
case. 



5 Appendix. Modified Novikov-Veselov equa- 
tions with periodic boundary conditions. 

In this Appendix we discuss the zero-energy spectral transform for the double- 
periodic Dirac operator (and the Generalized Weierstrass transform respec- 
tively) from the soliton point of view. This transform is naturally connected 
with a completely integrable hierarchy of integro-differential equations with 
2 spatial variables known as Modified Novikov-Veselov hierarchy (MNV). 

Formally, the results of this Section are not used in our proof of Theo- 
rem [I], but we hope they allow the reader to gain a better understanding of 
the problem. 

There is a rather big amount of papers dedicated to the periodic prob- 
lem for soliton equations (see for example the textbook pip . Nevertheless, 
the direct spectral transform for two-dimensional Dirac operator was never 
studied in such context in the literature known to us. Important properties 
of Bloch varieties for multidimensional Dirac operators were proved in JT3], 
but they are not sufficient for the purpos of integrating the periodic MNV 
equations. 

From the point of view of the Bloch theory the two-dimensional double- 
periodic Dirac operator is rather similar to the two-dimensional double- 
periodic Schrodinger operator. The fixed-energy direct spectral transform 
for last one was constructed by soliton methods by I. M. Krichever [|10J. This 
problem was studied more detail by J. Feldman, H. Knorrer and E. Trubowitz 
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see |3J). In Section |5.1| we describe the structure of the Dirac Bloch variety 



by methods analogous to |10[ . It is important to remark, that in spite of the 
similarity between these two problems we have to overcome some additional 
technical difficulties on this way. The asymptotic expansion of the Bloch 
variety gives us "quasi-local" MNV conservation laws. 

The zero-energy scattering problem for the two-dimensional Dirac oper- 
ator possesses an infinite-dimensional algebra of symmetries, generated by 
the MNV equations. They were constructed by L. V. Bogdanov in [0. In [0 
a generalization of the Miura transform was defined, and it was shown, that 
this transform maps the MNV equations to the Novikov-Veselov hierarchy 
(see Hi~9fl), associated with the fixed-energy two-dimensional Schrodinger op- 
erator. MNV equations in the space of functions, decaying at infinity were 
integrated by the so-called method of <9-problem in [fj] . Periodic MNV theory 
is discussed in Section FT3. 



In contrast with the one- dimensional soliton systems, the two-dimensional 
ones essentially depend on the boundary conditions. To define the periodic 
MNV hierarchy uniquely we have to fix some constants of integration. One 
of the simplest way to do it is to define the MNV hierarchy in terms of the 
so-called Cauchy-Baker-Akhiezer (CBA) kernel. This kernel was introduced 
for the Kadomtsev-Petviashvily hierarchy by A. Yu.Orlov and one of the 
authors (P.G.) in |J (see also 0). MNV hierarchy in terms of the CBA 
kernel is discussed in Section |572|. 

In this Section we always assume that U(z, z) is real, smooth and double- 
periodic (||). 

5.1 Asymptotic structure of the Bloch variety. 

For large Impi, Imp2 the structure of the Bloch variety can be studied by 
methods of perturbation theory. Following |l(J we start from the Dirac op- 



erator with zero potential U(z,z) = 0. The corresponding Bloch variety is 
a union of two Riemann spheres T^ = T\ U T^ , T[ = T 2 = CP 1 with a 
coordinate A and 

$A, z,z)=(l) e x ", A e rS 0) , $A, z,z)=(°\ e Xz , A G lf\ (41) 
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A pair Ai G Ij , A2 G T^ is called resonant if 

e Al - Az = 1, e Alf ~ A2T = 1, (42) 

and non-resonant otherwise. All resonant pairs are given by the following 
formulas: 

A ( m , n) = nmRer-Tm + .^ ^ = ^ m>n) g z _ (43) 

Imr 

Let us call a point A G I^ non-resonant if the pair Aer| and A G T 2 is 
non- resonant. The antiholomorphic involution cr# maps T 2 to I]_ thus it is 
sufficient to develop a perturbation theory only on T\' . 

Let s, R be some positive constants. Denote by 1^ R the domain, obtained 

from CP 1 by removing e neighbourhoods of all resonant points A^ and the 
disk |A| < R. 

Lemma 2 For any e > there exists a constant R(e) such that in the do- 
main T £ L, there exists an unique solution of the Dirac equation ([Jj with 
normalization ffity such that 

$(\ z+l,z + l) = e x+h ^(X, z, z), 
$(\,z + r,z + f) = e Xf+h ^ T tf(\, z,z), { ' 

where h(X) is uniquely defined under the condition that h(X) — > as A — > 00. 
The functions ip(X, z, z), and h(X) are holomorphic in A in the domain T £ L. 

A proof of this statement is to appear in a forthcomming paper by I. A. Taimanov 
and one of the authors (M.S.). It is rather long and rather technical. We do 
not want to present it in our text. 

Statement 2 The functions ip(X, z, z), h(X) defined in Lemma [| posses the 
following asymptotic expansions as A — ► 00 

' 1 1 <t>i{z,z) I <l>2(z,z) I <t>3(z,z) I <t>4,(z,z) I \ 

1 T i ~T \2 T \3 \4 1" . . . I 



X(\ 7 =\ _ p A(*-*i)+h(A)(*-*i) J- 1 A r ^?~ 1" ^3— -T — ^ 



A "+" A 2 ~+~ A 3 ~+~ A 4 ~+~ 



hi h z , h 5 



(45) 



fi/je Block variety V has a symmetry a : (wi,W2) — ► (wi , w^ )? °"(^) = ~~ A, 
t/itts a// even coefficients in fiffl) are identically zero.) 
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Unfortunately, at this moment we do not have any complete proof of this 
Statement. It is rather clear how to do it, but this proof needs a rather big 
piece of asymptotic analysis, and we are not in a position to do it presently. 
But we know, that it is fulfilled at least in two important specific situations: 

1. U(z,z) is an algebraic-geometrical (or, equivalently, finite-gap) poten- 
tial. It means, that the normalization of the zero-energy Bloch variety 
is algebraic (has finite genus). 

2. U(z, z) depend only on one real variable x — Re z. This fact was proved 



by one of the authors (M.S.) in fl6 |. Such potentials corresponds to 
surfaces of revolution. 

We would like to remark, that in |4j] a class of Riemann surfaces was intro- 
duced, which are in some sense similar to compact Riemann surfaces, and 
the zero-energy level of a two-dimensional Schrodinger operator belongs to 
this class. Therefore it is natural to expect that the zero-energy level of the 
two-dimensional Dirac operator also belongs to this class. 

To define the Modified Novikov-Veselov equations (MNV) and their con- 
servation laws it is sufficient to have a formal solution of the Dirac equation 
in the form ([45|), (|46|). Let us show that such solution exists and is unique if 
we assume that all (f>k(z, z), Xk(z, z) are bounded in the whole z-plane. 

Inserting ([45|) and 

MA) = y + X f + X f + ... (47) 

in (||) we get the following system of equations 

Xi(z,z) = -U(z,z), 
Xk(z, z) = -dzXk-i(z, z) - U(z, z)(j)k-i(z, z), k > 1, (48) 

d z (fr k (z, z) = U(z, z)xk{z, z)-h k - £*=i hj<\) k _j{z, z). 

We solve this system by induction. First we find Xi{ z i %)■> then (f>i(z, z), then 
X2{z, z), then (f> 2 (z, z) and so on. To find Xk( z , z) at each step we differentiate 
some double-periodic functions, obtained at previous steps. Thus they are 
defined uniquely and are automatically double-periodic. To find <fik(. z ,z) 
we have to invert the operator d z in the space of functions, bounded in 



the whole z-plane (any bounded solution <fik(z,z) is automatically double- 
periodic). This is possible if and only if the mean value of the right-hand 
side is equal to zero: 

fc-i 
« U(z, z)xk(z, z)-h k -J2 h j<Pk-j(z, z) »= 0. (49) 

i=i 

where: 

«F(z,z)»= f [ dtidt 2 F(t 1 + rt 2 ,ti + ft 2 ). (50) 

JO JO 

Thus at each step we find hk from ([19|) and then calculate <fik(z,z). The 
function <f>k(z,z) is determined by the system (|48| ) uniquely up to adding 
an arbitrary constant. This constant is fixed by the normalization condition 

Let us check that the function h(X) does not depend on the normalization 
point z\. If we change the point z%, then we change the integration constants. 
But these constants can be arbitrary shifted by multiplying the whole solution 
to a formal series in A 



4>l(z,z) I <t>2Jz,z) 



Xi(z,z) | X2Jz,z) , 



A A 2 



]_ _i_ <fti(z.*0 _i_ <h{z,z) _|_ 
A A 



(51) 

This multiplication does not affect h(X). 

We have proved, that the constants hi, h 3 , h 5 , . . .are completely deter- 
mined by the potential U(z, z) and all h 2 k = 0. Thus we have constructed 
an infinite sequence of functionals h 2 k+i[U]. The formulas for the first two 
of them are: 

h x = - «U 2 (z,z) » (52) 

h 3 = -« U{z, z)U- z - z (z, z) - (U 2 (z, z) + hx)V ls {z, z) », (53) 

where 

V 1 (z,z)=d;\U 2 (z,z)+h 1 ). (54) 

Let us point out that adding an arbitrary constant to Vi(z, z) does not affect 
h 3 . 

We have defined an infinite collection of functionals h 2 k+i[U], k — 0, 1, 2, . . 
The following Statement explains why these functionals are so important. 
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Statement 3 The quantities h^k+i [U] are conservation laws for the whole 
hierarchy of the Modified Novikov-Veselov equations (MNV). 

Assuming that Statement |^ is fulfilled we the prove this Statement at the 
end of the Section. 

It is well-known in the soliton theory that integrable systems with one 
spatial variable usually have infinitely many local conservation laws. For 
multidimensional soliton systems we normally have the opposite situation: 
almost all conservation laws are nonlocal. Let us briefly discuss the case of 
the MNV hierarchy. 

A functional Q[U] is called local if it possesses the following representa- 
tion: 

Q[U] =« q(U, U z , U z , U zz , U z - Z , U zz , ...)» (55) 

where the density q(. . .) depend only on U(z, z) and finite number of its 
derivatives. Of course hi[U] is local. The next conservation law h 3 [U] is 
nonlocal because the corresponding density depends on an auxiliary func- 
tion Viz(z,z), and to calculate Viz(z,z) we have to know U(z,z) on the 
whole z-plane. It it rather evident that higher functionals h2k+i[U] are also 
nonlocal. This nonlocality creates no serious problems if the potential is 
double-periodic, but it is very difficult to extend existing definitions to wider 
classes of boundary conditions. 

In [[H]] perturbation theory was developed also in neighbourhoods of res- 



onant pairs. It can be shown that for sufficiently large A the surface T is 
obtained from the T(0) by attaching small handles to the resonant pairs. We 
do not use this fact, thus we do not want to discuss it now. However we use 
the following property of T: 

Corollary 1 The Bloch variety F has two infinite points, corresponding to 
the points A = oo in T[ and T 2 ■ We denote them by oo + and oo_ respec- 
tively. 

This statement follows immediately from Lemma |2|. 

Further we shall use the following notation: A — ► oo + , where A is a point 
of T, always means that A tends to oo + in the domain F £ R , £ -,; the notation 

A — ► oo_ always means that A E T^ L. 
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5.2 Cauchy-Baker-Akhiezer kernel. 

To study symmetries of the soliton equation it is convenient to use the 
Cauchy-Baker-Akhiezer kernel (CBA) (see 0). In this Section we define the 
CBA kernel on the zero-energy Bloch variety of the two-dimensional Dirac 
operator. 

To start with, suppose (A, /j,) is a pair of points in T such that 

1. if)(i>,z,z) is non-singular at the points v = A and v = cr/i. (Let us 
recall, that outside v = oo the poles of ip(v, z, z) arose due to the 
normalization ( |T2| ) and do not depend on z and z.) 



2. At least one of the following conditions is fulfilled: 

Imp^A) — Imp!(/x) /0 or Imp 2 (A) — Imp 2 (/-0 7^ 0. (56) 

Then we may define u)(A, //, z, z) by: 

ui(\,(j,,z,z)= / du(X, fj,, z\ z'), (57) 

J oo 

where 

du)(X,fj,, z', z') = ^2 (A, 2 ; , z')ip2{crn, z', z')dz' — i>i(X, z', z')ipi(afi, z', z')dz' 

(58) 
The integral in (]57[) is taken along some path 7 in the ^-plane such that 



1. 7 connects the point z with 00. 

2. The form duj(X, /i, z', z') decays exponentially as z' — > 00 along 7. 

Condition ( p^ ) guarantees the existence of such path. 

From the Dirac equation (|4]) if follows that the form duj(X, fi, z',z') is 
closed that the integral (0) is well-defined. 

The next step is to show, that for a fixed \x our function u)(X,[j,,z, z) is 
meromorphic in A on V outside the points oo + and oo_. 

Suppose, that a pair (A, //) satisfies the following condition, which is much 
weaker than (156T). 



2'. At least one of the following combinations 

— (pi(A) - PiQjt)) or -— — - (p 2 (A) - P2(A*)) is non-integer, (59) 

Zn Zti\t\ 
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or, equivalently, at least one of the following inequalities is fulfilled 



Wi(X)w 1 l {n) ^ 1, or w 2 (X)w 2 1 (fi) ^ 1. 
Assume that Wi(\)w~ (fi) ^ 1. Then we have either 

o)(A, fj>,z,z) = / uj x (X, n, z + t, z + t)dt 



or 



w(A,/i,z,z) 



cu x (A, n,z + t,z + t)dt, 



(60) 

(61) 
(62) 



where 



CJ X (X, fj,, z', z) = if> 2 (\, z', z)%l) 2 {a^ z', z) - ^i(A, z', z^^a/j,, z', z'), (63) 



and t ER. 

If Wi(A)wf 1 (/i) > 1 from (|6~lD we get: 

w(A,£t,3,*) = 
= f / + 7 + 7 + • • • J W*(A, M, 2 + t, 2 + t)dt = 



(64) 



mJB) 



Wl(A)-1Ul(/i) 



^y / W x (\, jJ,,Z + t,Z + t)dt 



Similarly if W\{\)w 1 1 ( / u) < 1 from (|62]) we get: 

w(A, /i, z, 2) = Ml( ff l(/i) / £*(A, fi,z + t,z + t)dt 



(65) 



The formulas (p4|) and ( p5|) define meromorphic continuations of the function 
u)(A, /i, z, z) to the whole surface Y. To conclude the proof, it remains to note 
that formulas ( |64|) and (|65|) coincide. 

If W2(A)w 2 ^ 1 (/i) 7^ 1 then the integration path can be chosen along the 
line z' = z + rt, and we obtain a new representation for the same function 



Cj(X, fi, z, z) = , llo( ^L. A I w 2 (A, n,z+rt,z + Tt)dt, 



W2(X)—W2(n) 



(66) 
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where 
uj 2 (X, n, z' , z') = Tip 2 (X, z' , z')ip 2 (an, z' ', z) - fijJi(\,z,z)ilJi(afj,, z',z). (67) 
Now we are in position to define the Cauchy-Baker-Akhiezer kernel 

w(A, fi, z, z) = -^-——, zr — w(A, H, z, z) (68) 

2n< u x ([i,[i,z,z) > x 

where 

< uj x (fi, n, z, z) > x = / Cu x (\, fi, z + t, z + t)dt. (69) 

Jo 

It is a simple exercise to check that the function < uj x (fi, fi, z, z) > x does not 
depend on z and is even in /i, i.e. 

< Q x (fx, fi, z, z) > x =< uJ x (a/i, afi, z, z) > x . (70) 

c (±) c (±) 
< w x (/i, //, z, z) > x = =Fl + -^- + -^j- + . . . , /i -»• oo± (71) 

Lemma 3 The Cauchy-Baker-Akhiezer kernel cu(X, n, z, z) defined above has 
the following properties: 

1. For any fixed z, oo(X, /j,,z, z) is a meromorphic function of X and a 
meromorphic 1-form in /j, (both on the finite part ofT). 

2. For generic ji, cu(X, /i, z, z) has poles at the poles ofip(X, z, z) and at the 
point X = fi. It is holomorphic outside these points at the finite part of 

r. 

3. For a generic X and \x — > A 

u>(X, \x, z, z) = + regular terms, (72) 

2m n — X 

or, equivalently, 

u(X,ji,z, z) = 1, (73) 



where S is a small contour, surrounding the point X. 
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4- Let \x — > oo + and let X be a fixed point in T. Then we have the following 
formal expansion: 



Uj(X,Ll,Z,z) 



where 



dfj, 

2ixi 



M\z,z) + ^B^[d 2 ]M\z,z) 



/' 



k=2 



fe-2 



II* 



-fiz—h(fi)z 



(74) 
(75) 



1=0 

are differential operators in z of order k — 1. The coefficients v kl (z, z) 
do not depend on X and ll, are differential polynomials of the asymptotic 
expansion coefficients <j)j(z, z), Xj{ z i z)> j < k. Each function v kl (z, z) 
depends also on a finite number of constants h 2 j+i, cL . 

Similarly for il — > oo_ ; we have: 






^(a,v-) + j4; ) [#(a 1 ^i 



// 



k=2 



where 



k-2 



R k - ) m = d k z -^ + Y J v k - l \zrz)d l z . 



ti« 



(-)/.» =^ 



-/iz—h(fi)z 



(76) 
(77) 



1=0 



u(X, li, z + 1, z + 1) = w 1 (X)w 1 1 (/i)u;(A, fj,, z, z), 
w(\, Li, z + r, z + r) = w 2 (X)w 2 ~ 1 (li)uj(X 1 li, z, z). 



(78) 



Statement 4 If Statement [| is fulfilled, then the CBA kernel has the fol- 
lowing additional properties: 



1. The formal expansions (\74i), ( \76j) are asymptotic. 

2. Let /^(A, z, z) be the following formal series in X 



/™(a,m) = < £ i£^l\e^ 

k=-N 



x k 



(79) 
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Then 



2iri res 



u(X,iJ,,z } z)f {+) (iJ,z,z) 



J2 tPj0 i c ( l 



fl=OC^ 



k=-N 



x^ 



„Az 



as X — > oo + and 

2ni res 



u,{\,IA,z,z)f™(iM,z,z) = o(\)e»*, 



(80) 
(81) 



as A — ► oo_. 

Similarly if f(~'(\, z, z) is the following formal series in X 

f^(X,zrz) = lf £^0\e» 



k=-N 



x k 



(82) 



then 
2iri res 



U,{\,p,Z,z)fM(vL,Z,z) = { £ ^ ! ( f ' ^ + Q ( 1 



H=oo- 



k=-N 



A* 



A 



_A.~ 



as A — > oo_ ana" 



27ri res 



u(X,n,z,z)f { \fi,z,z) 



O i j) e As , 



(83) 
(84) 



as X — ► oo_ 



Remark 4 Tae formulas (j^W, (ffl\), (BS), ( ^i ) require some comments. On 
the left-hand side we have formal series in \x, thus the analytic definition of 
the residue as an integral does not work. Fortunately the terms, containing 
exponents in fi annihilate each other and we can define the residue as the 
coefficient ofl/ji: 



res 



-j-dfi = at 

^=°° k=-N A* 



(85) 



Proof of Lemma 0. 
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1. We have constructed u)(A, //, z, z) as a meromorphic function. The func- 
tion < uj x (lj,, (jl, z, z) > x is meromorphic in li on the whole T and dpi(fi) 
is a meromorphic 1-form on T. Thus (65) gives us a meromorphic 



function with appropriate tensor properties. 
2. From Lemma |2| it follows that for generic ii 

(a) ^(i/, z i z) is nonsingular at the points v = /i, v = aii, 

(b) < u} x (fi,fj,,z,z) > x ^ 

(c) wi(A)wf (n) = u7 2 (A)w^ 1 (/i) = 1 iff A = li. 

If fi fulfill these conditions, if A ^ ii and if ip(X, z, z) is nonsingular, 
then formulas (|65|), (|66|), (|68|), define a nonsingular function of z, z. 



3. Let A — -> /i. To calculate the asymptotics of the CBA kernel we expand 
the function W\(X) in (|>|) in a series in A — fi using ([IB]), (|i4|). We see, 
that we have a first-order pole and the normalization in (^) is chosen 
so that the residue is exactly 1. 

4. To calculate the asymptotics of u(X, li,z,z) as /i — * oo we substitute 
the asymptotic expansion (|45| ) to (|57|) and use the following asymptotic 
formula: 

Azjva „/ ='\J J , „Af /\ / -A j-/ Air^/ i\/t-1^2 #(A, £,^) 



| e A V(A, z', z')^' + e Af ^(A, z', z')dz' = e A " E(" X 



A fc ' 

(86) 
where the functions f(X,z',z / ), g(X,z',z') are some asymptotic series 
in A 

/(*,*,*) = !, — TT-, 9{\z,z) = 2^ — -— (87) 

fc=0 A fc=0 A 

such that 

3, (e A 7(A, z', z')) = d z (e x "g(X, z', z')) , (88) 

and all expansion coefficients fk{z, z), gk{z, z) are smooth functions, all 
derivatives of these functions are bounded in the z-plane. 

The proof of the formula (|86|) is rather standard and we do not repro- 
duce it here. 
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5. The last statement of the Lemma follows directly from the transforma- 
tion rules for the Bloch functions. 



Proof of Statement |. Assume for defmiteness, that A — ► oo + . If \i 
then: 

1 ^4l + \z,z) 



From (jS9| ) it follows, that 






AV 



,(A-/i)z 



dfi. 



2iri res 



^(A,/i, z, z)f i+ \fi,z,z) 

fl=OD + 



res 



dfi 



,{X-ti)z f(+)i 



/i=0O-|- fJ' A 

Applying the well-known formula: 



r .;~fW(( iiZl z) + 0\-\( 



A 



a; 



res 



/i=oo 



dfi 
fi — X 



/'' 



X n n > 
n <0 



we get the right-hand side of (j80|) . 



If fi — > oo_, then using (|7^) we get: 



ij>0 



wl +_) U,z) 



AV 



\Z — [LZ 



dfi. 



Formula fl84|) follows automatically from from (p2|) . 

Formulas (B3f), (BT) are proved exactly in the same way. 



oo_ 



(89) 



(90) 



(91) 



(92) 



Remark 5 A formula, representing the Cauchy kernel on a Riemann surface 
as a semi-infinite sum of quadratic combinations of eigenf unctions first arose 
in the article / [77| / by I. M. Krichever and S. P. Novikov. In j71|/ the spatial 
variable x was discrete. In j^J the spatial variable was continuous, and the 



CBA kernel was defined as an integral similar to ffi 7 /). 



Remark 6 A form, similar to u(X, fi,z,z), but with integration path, start- 
ing from a finite point Z in the z-plane, arises in the theory of finite Darboux 
transformations (see ]73|/, formula (6.1.18)). In fact, the same object arose 
in the Generalized Weierstrass map. We do not discuss this analogy in our 
text, but it seems possible that this analogy has some deep roots. 
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5.3 Modified Novikov-Veselov equations. 

Lemma 4 Let n be a generic point in Y . Define a deformation <jw of the 
function i[>(\, z, z) associated with this point by: 

<5 (/i Vi(A, z, z) = w(A, fi, z, ztyxtji, z, z) + a(X, n)i[>i(\, z, z) ,_ 

S^ip 2 (X, z, z) = Q(\, h, z, z)ip 2 (fi, z, z) + a(X, ^^(A, z, z), 

where a(X, (x) is a meromorphic function in X, uniquely fixed by the following 
requirement 

c^Vi(A, z, z) + 6^ 2 {\ z, z)\ = 0. (94) 

Then (ffi\) is a deformation of the Bloch function, corresponding to the fol- 
lowing deformation of the Dirac operator 



5 {fl) L 
where 



-S[^U(z,z) 

S^U(z,z) 



(95) 



6f ) U(z, z) = ^i(/i, z, z)tp 2 {o-fi, z, z) , gg , 

S^U(z, z) = V> 2 (/i, z, z^io-fi, z,z). 

For generic fi these deformations result in non- self- adjoint Dirac opera- 
tors (8\ U(z,z) 7^ 82 U(z, z)) with complex-valued potentials . But the the 
Bloch function and the Bloch variety are well defined for such Dirac opera- 
tors. 

Proof of Lemma |j. A simple direct calculation shows, that 

(L + eS^L) (tf(\, z, z) + e6Mtf(\, z, z)) = 0(e 2 ), (97) 



Thus (|93|) defines deformations of the eigenf unctions, corresponding to (|9 



From (ff8|) it follows that the new eigenf unctions satisfy (11) with the same 
multipliers W\(X), w 2 (X) as the old ones. Thus the new eigenf unctions are 
defined on the same curve T and have the same periodicity properties. The 
last property plays the key role, when we prove below that the functionals 
h 2 k+i[U] are invariant under some deformations. 

The deformations of the Dirac operator, generated by all 5^' form a linear 
space. Deformations preserving the class of self-adjoint Dirac operators with 
real potentials are the most interesting ones. Let us check that the subspace 
of such deformations is sufficiently large. 
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Lemma 5 Denote by A^> the following linear combination of deformations 
SO*).- 

A^ = — + — — - (98) 

< uj x (fi,fi, z, z) > x <u x (8fi,6fi,z,z) > x 

(recall that < o) x (/i, n,z,z) > x is an even function in \x (see ffi\)) and does 
not depend on z, z.) 

Then A^) acts on the space of self-adjoint Dirac operators with real po- 
tentials, i.e. 



A^ ] U(z, z) = A^ ] U(z, z) = A { ^U{z, z) = A^U{z, z) 

(99) 

and the Bloch variety T[U] defined in Section [3| is invariant under these 
deformations. 



Tpi(lJ,,z,z)ip2(crfi,z,z)+ipi(<jfi,z,z)4>2(lJ-,z,z) | i/ji(9fi,z,z)4>2(c8lJ;Z,z) J t-Tpi((T8fi,z,z)ip2(8fi,z,z) 

<UJ x (fl,fJ,,Z,z)> x <LU x (8fl,9lI,Z,z)> x 



Proof of Lemma ^. Formula ( |9"9"D follows directly from fl9"E| ) and (fJSI). The 
Bloch function ijj(\,z,z) for real U(z,z) has symmetry property fllTD, thus 
the right-hand side of (^Uj) is real. As we pointed out in the previous Lemma, 
all deformations generated by S^' do not change the Bloch variety T[U] and 
the multipliers w%, w%. This completes the proof. If Statement ^| is fulfilled, 
then these flows do not change h 2 k+i[U}. 

Equations fl99|) are essentially nonlocal and rather complicated. Namely 
the right-hand side is expressed in terms of the Bloch function, and it is 
difficult to calculate Bloch solutions either analytically or numerically. For- 
tunately the space of deformations generated by all A^ contains simpler 
equations such that the right-hand side can be expressed via U(z, z) in terms 
of quadratures. 

Let n — > oo. If Statement is fulfilled we may expand A^ to the 
following asymptotic series 



(We write the term d/j,/dp(fj.) to gain some standard normalization of the 
MNV hierarchy. If we omit this multiplier, our expansion coefficients will be 
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linear combinations of Novikov-Veselov generators with constant coefficients, 
which in most situations is not essential). 

Any K 2k+1 [U] is a quadratic polynomial of (f>i(z,z), xi{z,z), I = 1, . . . , 2k, 
with coefficients depending on the h 2 i-i and c 2 \, I = 1, . . . , k, where the c 2 ; 
are coefficients of the asymptotic expansion ([n]) . 

For any odd integer 2k + 1, k > we have the following pair of flows on 
the space of real double-periodic functions: 

dU{zrzMk+l) - = 2ReK 2k+1 [U]- (101) 



dt 2 k+\ 
dU(z, z, i 2k+1 



dt 



2k+l 



2lmK 2k+1 [U}. (102) 



Definition 1 The equations t \101\ ), ( \103j ) are called the Modified Novikov- 
Veselov equations (MNV). 



Remark 7 To define the MNV flows associated to ( \101\ ), (10i) it is suf- 



ficient to have a formal expansion for the Bloch function and the function 
h(X). Thus these flows are well-defined for any smooth potential. But without 
Statement [| we can not prove that they conserve the functional h 2k+ i[U}. 

Using this definition of the MNV hierarchy we may immediately prove 
Statement || All functionals h 2k +i [U] defined in Section |3| are integrals of 
motion for all flows A^, thus they are conservation laws for their expansion 
coefficients. 

The representation ( p.01| ) may look rather unusual. To check that our 
definition of the MNV equations coincides with the standard one, let us 
calculate the deformation of the Bloch function, corresponding to the flows 



( 101 ). To gain a standard answer we shall use a normalization of the Bloch 
function different from th eone used earlier. Instead of (|94j) we assume that 
the function a(A,/i) in (|93|) is identically zero. 

Statement 5 Let the Statement || be valid. Then the deformation of the 
Bloch function corresponding to the flows ( \101\) has the following represen- 
tations: 

dip(\,z,z,t 2k+1 ) _ 

dt 2k+ \ 
30 



2ni < res 



+ res 

/i=00-|- 



U}(X,fi, Z, Z,t 2 fc+l,)/i ' + ip(pL, Z,Z,t 2 k+l), (103) 



/j=oo_ 



A 2fc+ V(A,^,W) + 



O Qj e A * as A -»• oo_ 

O h) e Xz as A ->• oo. 



(104) 



= {df +1 + Sf +1 + E Wifo *)# + Wfo ^)4} #A, z, z, t 2k+1 ), (105) 
Proof of Statement [5], given the validity of Statement |2|. By definition 



idp(n) 



A (/ V(A, *, z) 



= m 
+7ri 



u(X, \X, z, z)ijj(fi, z, z) + \u{X, — //, z, z)ij)(—fji, z, z) 
u(X, ft, z, z)if)(ft, z, z) + [oj(X, —ft, z, z)ij)(—ft, z, z) 



+ 



p— >oo+ 



(106) 



(in this formula we use o\i = —/a, 9fi = —ft). 

The residues in (|103| ) are exactly the coefficients of the terms dfi/n 2k+2 
and dft/ft 2k+2 respectively. Thus these coefficients gives us the action of 
K 2 k+i[U] and K 2k+ i[U] on the Bloch function. 

Formula ( |104j) follows directly from ( |103| ) and formulas ( J30|) - (|8~4T ). 

To prove ( |105| ) let us substitute the asymptotic expansions (|74]) , ( 
QID3D . We get: 

dipi(X,z,z,t 2k+1 ) 



111 



dt 



2k+l 



>(+) 



2A- 



(+)[ 



R 2 y + i%\ + e <hk+x-i{z,z)Rr[d*] \M^z,z,t 2k+1 )+ 



2k 



,(") 



?(-) 



(107) 



+ E X 2k ii-i(z, z)R\ > [d z ] ^ 2 (A, z, z, t 2k+1 ) 



.1=0 



= dl k+l il)i(X,z,z,t 2k+ i) + U(z, z)d 2 ip 2 (X, z,z,t 2k+1 ) + lower order terms = 

(108) 
= [df +1 + d 2k+1 } Vi (A, 2, z, £ 2fc+1 ) + lower order terms, (109) 

(the x\ \ \Z, z) denote the expansion coefficients of the function ipi(X, z, z) at 
the point A = oo_). A similar calculation shows the following relation: 

d^(A, z, z, t 2fc+1 _) = | Q2fe+1 + a2fe+1 | ^ ^ _^ t ^ +i) + bwer order termg _ 

(110) 



dt 



2k+l 
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Statement [5] is proved. 

From (|105|) it follows that the MNV equations are the compatibility con- 
ditions for a pair of differential operators on the space of zero eigenfunctions 
of L. It is well-known in soliton theory that such compatibility conditions 
are equivalent to the existence of standard L — A — B representations (see 

!)• " 

Remark 8 In this article we constructed a soliton hierarchy in terms of the 
Cauchy-Baker-Akhiezer kernel. 

Using the CBA kernel we construct, in fact, a much wider hierarchy, in- 
cluding essentially nonlocal equations. All these equations preserve the spec- 
tral curve. In Section we show that in fact the deformations corresponding 
to conformal transformations of the Euclidean space, lie in this wider hier- 
archy. 

In terms of U(z, z) this hierarchy looks rather unnatural. But we may 
treat it as a system of differential equations on a bigger collection of functions 
simultaneously; namely we may consider the potential U(z, z) and the wave 
function in a finite number of fixed points on the spectral curve as unknown 
functions, connected by the Dirac equations. Similar systems associated with 
one- dimensional soliton equations were discussed in the literature (see jT^ / 
and references therein) from both a mathematical and a physical point of view. 
In ]7o|/ it was shown, that starting from the KP equation we get a hierarchy 
naturally containing many other well-known soliton systems. 
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